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Abstract 

In this paper we investigate predictability of electricity prices in the Canadian provinces of Alberta and Ontario, 
as well as in the US Mid-C market. Using scale-dependent detrended fluctuation analysis, spectral analysis, and the 
probability distribution analysis we show that the studied markets exhibit strongly anti-persistent properties suggesting 
that their dynamics can be predicted based on historic price records across the range of time scales from one hour to 
one month. For both Canadian markets, the price movements reveal three types of correlated behavior which can 
be used for forecasting. The discovered scenarios remain the same on different time scales up to one month as well 
as for on- and off- peak electricity data. These scenarios represent sharp increases of prices and are not present in 
the Mid-C market due to its lower volatility. We argue that extreme price movements in this market should follow 
the same tendency as the more volatile Canadian markets. The estimated values of the Pareto indices suggest that 
the prediction of these events can be statistically stable. The results obtained provide new relevant information for 
managing financial risks associated with the dynamics of electricity derivatives over time frame exceeding one day. 
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1. Introduction 

The modern electricity market is not only a system for arranging the purchase and sale of electricity using supply 
and demand to set the price, but, for most major grids, is a basis for electricity derivatives, such as electricity futures 
and options, which are actively traded. The practical significance of this part of the market is increasing as is the 
importance of the related scientific research [1, 2, 3]. The markets of electricity derivatives have developed as a result 
of the liberalization and deregulation of electric power systems around the world. Deregulation, introduced initially to 
reduce and simplify the control of the business in this field, had a final goal to reach financial efficiency of electricity 
markets [4, 5]. However, electricity is unique as it is a non-storable commodity, and the markets remain extremely 
inefficient [6, 7]. 

Electricity prices are not a result of long-term but instant, usually on an hourly interval, balance of supply and 
demand. Moreover, as a consequence of the complexity of a wholesale electricity market, it can show an extremely 
high price volatility at times of peak demand and supply shortages. This price spikes are hard to predict and financial 
risk management is still a high priority for participants in deregulated electricity markets due to the substantial price 
and volume risks that the markets can exhibit [8, 9, 10]. 

The problem of predictability of electricity prices in deregulated markets has been considered in many previous 
studies (for instance, [11, 12, 13, 14]). The values of prices can vary by a factor of 100 over a time scale of just several 
hours. These dramatic changes tend to occur in a seemingly spontaneous fashion which is sometimes erroneously 
interpreted as a signature of a random uncorrelated process (see for example [15]). A more detailed mathematical 
analysis reveals nontrivial auto-correlations in these sudden price jumps [16, 17, 18, 19] which indicate a possibility 
of prediction of electricity price movements based on the information on their historic evolution [7]. However, it is a 
widely recognized fact that price fluctuations in energy markets display heavy distribution tails [16, 17, 20, 21, 22, 23] 
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causing substantial difficulties in building quantitative forecasting models of price behavior. Less attention has been 
paid to the analysis of temporal patterns underlying the observed statistical structure of electricity markets and as a 
result modeling of their dynamics is still in its infancy and is typically limited to day-ahead models [2, 24, 25, 26, 27]. 

In this study, we take a few next steps toward answering fundamental questions related to the predictability of 
electricity prices. First of all, can deregulated electricity markets reach the state of efficiency with the Hurst index 
value close to other well-known markets, or this state is not reachable in the usual sense [28, 29, 30]? This is crucially 
important because if the electricity markets are inherently inefficient, the forecast can be built at various time scales. 
In this context, the inefficiency means that price history is relevant to the future price changes and can be used for their 
forecasting [31]. The problem with Pareto-type statistics plays a special part, because not all heavy distribution tails 
can be approximated by a single probabilistic model. They can include several dynamic ranges described by distinct 
Pareto exponents. If such markets are predictable in principle, there might be particular price intervals for which the 
forecast is statistically stable, and these intervals are important to identify. 

Another central question related to the predictability of electricity prices is how universal can be a model of 
electricity price behavior across different markets. In the present work, this question is addressed in frames of a 
quantitative analysis of electricity prices in three independent markets with different levels of liberalization - Alberta, 
Ontario (Canada) and Mid-C (USA) markets. 

Dynamical and statistical properties of price fluctuations are investigated using several methods. First, we evaluate 
correlations in price dynamics across different time scales using the method of scale-dependent fractal exponent [7] 
obtained from detrended fluctuation analysis (DFA) [32, 33, 34, 35], We also use the Fourier spectral analysis to 
identify cyclic components in the electricity price dynamics, as well as the Pareto probability distribution analysis for 
testing the stability of statistical moments of the studied data. Spectral and DFA analysis results show no evidence 
of informational efficiency of electricity price fluctuations at any time scale. All three markets demonstrate different 
levels of inefficiency which could reflect their different sizes and structural diversification. Price movements in these 
markets are strongly anti-persistent [28], Together with Pareto analysis results, this anti-persistence indicates that 
electricity price movements can be predicted based on historic price records. 

Next, we verify the possibility of price forecasting using phase diagrams representing the correlation of previous 
and current price increments. According to our results, the diagrams have a complex asymmetric shape revealing 
three basic scenarios of price movements. These scenarios remain the same for price movements at different time 
scales, from one hour up to one month, and are found to represent strongly volatile market conditions. Based on these 
results, we show that price fluctuations in deregulated electricity markets are predictable by their nature. Our findings 
lay a foundation for future mathematical description of multiscale dynamics in deregulated electricity markets. 

The plan of the paper is the following. Section 2 contains a detailed description of the analyzed data sets. Section 
3 describes main results of our statistical analysis demonstrating the possibility of electricity price forecasting. This 
possibility is explored further in Sections 4. Section 5 provides a brief summary of our study. 

2. Data 

As an outcome of the liberalization policies pursued in several countries from the 80s on, the so called day-ahead 
electricity market provides economists with a very challenging phenomenon. Electricity cannot be economically 
stored, which implies that demand and supply must be continuously balanced, so that the market price mainly reflects 
the demand and supply conditions prevailing at the very moment it has to be delivered to final users. Then, rather 
complex market systems have been set up, with the aim of reaching a reasonable trade-off between economic efficiency 
and system reliability. These systems are built around a market operator, whose task is to manage uniform-price, 
sealed-bid, bilateral auctions in order to construct aggregate demand and supply curves, and to determine equilibrium 
prices and quantities. The knife-edge character of such a price setting mechanism is fatherly pushed to the extreme 
by a very low price elasticity of demand, and by technical constraints which time by time lead to network congestion 
(see e.g. [23, 36] and references therein). 

The data studied in this paper consists of hourly real time pool electricity prices in Alberta, posted by the Alberta 
Electric System Operator (AESO), and Ontario, posted by the Independent Electricity System Operator (IESO). The 
data cover the period from May 1, 2002 to June 6, 2009. In addition to these Canadian markets, the Mid Columbia 
(Mid-C) market has been considered during the time interval 1 July, 2001 to 31 Oct, 2006. For each of the three 
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Figure 1: Time series of hourly electricity prices in Alberta (left), Ontario (center) and Mid-C (right) markets. From top to bottom: all hourly 
prices, on-peak prices, and off-peak prices. Alberta electricity prices demonstrate significantly higher fluctuations than those in Ontario market, 
plotted on the same vertical scale. Fluctuations of electricity prices in Mid-C have twice as low amplitude as that in Ontario, and about 5 times 
smaller than in Alberta. 


hourly data sets, two secondary time series consisting of electricity prices during on- and off- peak hours have also 
been examined. Figure 1 shows the time series under study, including the original data and their on- and off- subsets. 
Note that all plots contain numerous spikes with irregular timing and amplitude. 

Alberta and Ontario are the only two Canadian provinces where wholesale electricity markets are fully deregulated 
[13, 37], Alberta’s market is dominated by fossil fuel generation and as such follows more closely the price of natural 
gas. Ontario’s generation involves about 50% of nuclear and 25% of hydro power [38, 39] enabling a more stable 
price behavior [40], The average level of volatility of electricity prices in Alberta is about twice as high as in the 
Ontario market. 

The Mid Columbia electricity market is not as deregulated as Alberta and Ontario are [41]. It is not a centralized 
power market, but it is a trading hub where power is bilaterally traded among utilities and marketers. The Mid-C price 
hub is a reference price for the Pacific NW region, which consists of Washington, Idaho, and Oregon. In this region, 
large utilities own generation and serve load under regulated rates. The generation is primarily hydro and the region 
typically exports to British Columbia and California [7, 42]. For these reasons, Mid-C prices are significantly less 
volatile than those in either Canadian market. 

3. Statistical signatures and predictability 

3.1. DFA and spectral signatures 

For testing the informational efficiency of electricity price fluctuations, multiscale correlations of price dynamics 
were evaluated across different time scales. Two complementary approaches were used to achieve this goal - the 
scale-dependent DFA and the Fourier spectral analysis. 

The former of the two approaches has been first introduced in [7], In contrast to previous methods manipulating 
with average scaling exponents characterizing broad scaling ranges, we investigated the distribution of local DFA 
exponents over all time scales involved. This approach was shown to be the only suitable when the signals under 
study encompass qualitatively different types of behavior including random price movements, cycles, and spikes. 
Using the DFA as the base algorithm is justified by the presence of multiscale trends in the electricity data [29, 43]. 
The scale-dependent version of this algorithm presented in [7] enables the investigation of complex types of nonlinear 
behavior of financial and economic indicators by providing detailed information on the distribution of correlations 
over different scales, and is especially useful for quantitative analysis of market efficiency. 
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Figure 2: Dependence of the detrended variation F (eq. (2))and the local scale-dependent DFA slope a on the time scale n for all hourly, on- and 
off-peak electricity prices in Alberta (left), Ontario (center), and Mid-C (right) markets. The presented statistics reveal complex correlated structure 
of price movements with quasi-periodic components associated with daily and weekly cycles. In all presented data sets the scale-dependent DFA 
exponent is significantly below the level 1.5 defining the state of informational efficiency, which provides an opportunity of forecasting the prices 
over wide ranges of time scales. 


The DFA technique was applied to the time-integrated signal 

k 

y(k) = ^ (x(t) - (X)) , (1) 

t= 1 

in which (x) is the average value of the hourly electricity price x and k - 1, ...,1V, where N is the number of points in 
the studied time series. The integrated signal y(k ) was devided into M = N/n non-overlapping subintervals of equal 
length n ranging between 4 and 720 hours. The boxes were indexed by m = 1,..., M and their starting times were 
labeled by k nm . For every box, the least square regression line y„ m (k) representing the local linear trend in that box was 
fit to the data. Using these fits, the integrated series y(k) was locally detrended and the root mean square fluctuation 
of the resulting detrended signal was calculated. The described calculation was repeated for each of the M boxes and 
the resulting values were averaged to obtain the characteristic dependence on the time scale: 


F(n) = 



( 2 ) 


For a fractal (self-similar) financial time series x(t), the power-law relation between the root mean square fluctua¬ 
tion F and the time scale n is expected [44, 43]: 


F{n) ~ n a , 
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in which a is the DFA scaling exponent [32, 33], 

Note that our definition of the DFA exponent differs from that used in the majority of other studies of price 
fluctuations operating with logarithmic price returns, in which case a serves as a proxy to the Hurst index H. Because 
the time series of the electricity prices are quite spiky (especially in Alberta), the DFA exponents evaluated using the 
standard return-based approach would be close to zero, making quantitative cross-market comparisons statistically 
unreliable. To improve the accuracy of our analysis, we apply the DFA method directly to the time series of the 
electricity prices rather than to the price returns, which is consistent with the original formulation of the method [33]. 
Since we omit the step of calculating price returns, the resulting DFA exponent is greater than that derived from the 
returns by one, and so a — H + 1. 

Keeping this relationship in mind, efficient financial markets described by random Brownian walk models [45, 29] 
exhibit the values of a which are close to 1.5 signaling the absence of correlations between the price increments [28]. 
It has been shown [7, 46] that deregulated electricity markets do not satisfy this condition and in that sense are not 
efficient. In addition, a values of electricity markets tend to vary with scale reflecting complex structure of the price 
dynamics involving quasi-periodic cycles and random variations [7]. 

Fig. 2 shows scale-dependent behavior of the DFA functions F(/i) characterizing studied electricity markets. The 
local values of a exponents were computed within narrow n intervals of exponentially increasing width ensuring a 
uniform binning on a logarithmic scale. It can be seen that Canadian markets obey the anti-persistent condition a < 1.5 
across all temporal scales. The Mid-C market shows a persistent behavior characterized by a > 1.5 at n < 12. This 
indicates that the daily cycle in this market dominates random fluctuations leading to statistically significant trends 
over time spans shorter than the half-day interval. 

All three markets exhibit anti-persistent regimes at n > 24 and operate in a strongly inefficient state, confirming 
our previous result obtained for daily electricity prices [7]. Table 1 provides average and minimum DFA index values 
for the studied data, along with spectral analysis results discussed below. 

It can be noticed that the drop of a associated with the daily periodicity is shifted toward larger relative to n = 24. 
This shift may be caused by an interplay between the saturation of F{n) and a decrease of the trends at time scales 
greater that the cycle period. Qualitatively similar, albeit less pronounced shifted signature is observed near the weekly 
periodicity. 

The Fourier spectral analysis was used for a more accurate description of periodic components of price dynamics, 
and also as an independent test for its informational inefficiency. The power spectrum S (/) is obtained from the 
Fourier transform Ay of the price signal x(t) defined in the continuous limit as 

Xr(f) = f x(t)e- a « f, dt, S{f)=l-X t X* t , (4) 

Jo 1 

where Xj is the complex conjugate of Xj and T » 1 // is the length of the time interval of the analysis. The method 
was implemented using the standard fast Fourier transform algorithm. For a self-similar signal the spectral power 
scales with the frequency f as S (/) ~ f where the index f5 is related to the DFA index through/? = 2ar — 1, so that 
/? = 2 is the case of the efficient market (Fig. 3). 

The obtained Fourier power spectra are consistent with results of the DFA discussed above (Fig. 3), and reveal a 
strongly inefficient behavior of the electricity markets. In addition, spectra analysis demonstrates a relatively stable 
daily cycle as well as a less stable but evident weekly periodicity. The average f5 estimates are statistically significantly 
below the value 2. The values a t heor calculated based on the spectral exponents are in an approximate agreement with 
the average a values obtained directly from the DFA (the first and the last columns of Table 1). 

Both methods clearly show that there are strong anti-persistent correlations between the values of electricity prices 
at all time scales for all-hour time data as well as for the subsets of the on- and off-peak data. This is evident from 
the analysis of the maximum DFA index values (second column it Table 1) representing the largest a across the entire 
range of n scales for each market. In most of the data sets, [a(n)] max <1.5 meaning that at no scale the markets 
become efficient or persistent. The only exception from this tendency is the Mid-C market at n < 12 hours showing 
greater than 1.5 DFA index values for both all-hour and on-peak prices. This effect is absent in the off-peak Mid-C 
data suggesting that the short-term informational persistency of this market is a footprint of the high demand periods. 

Despite the significant variability of the DFA exponents across temporal scales, the average a values of all-hour 
prices in Alberta and Mid-C markets shown in Table 1 are in an agreement with those describing daily prices in these 
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Figure 3: Fourier power spectra S(f ) for all-hour, on- and off-peak electricity prices in Alberta (left), Ontario (center), and Mid-C (right) markets. 
In all data sets the average spectral exponent ft is lower than 2, which confirms the conclusion of the DFA analysis regarding the inefficient behavior 
of the studied price dynamics. Canadian markets are characterized by significantly smaller values of ft compared to the American Mid-C pool, and 
therefore have a higher degree of inefficiency. 


6 




























Table 1: Average values of DFA and spectral exponents of hourly electricity prices. 


Time series 

< a(n) > 


P 

& theor ~ 

08+0/2 

Alberta, all-hours prices 

0.87 + 0.13 

1.18 + 0.02 

0.93 ± 0.07 

0.97 

±0.54 

Alberta, on-peak prices 

0.87 ± 0.08 

1.09 + 0.01 

0.94 ± 0.06 

0.98 

±0.53 

Alberta, off-peak prices 

0.83 ± 0.09 

1.10 ±0.04 

0.81 ± 0.07 

0.91 

±0.54 

Ontario, all-hour prices 

0.87 + 0.13 

1.16 ±0.01 

0.92 ± 0.06 

0.96 

±0.53 

Ontario, on-peak prices 

0.86 ± 0.06 

1.03 + 0.04 

0.83 ± 0.06 

0.92 

±0.53 

Ontario, off-peak prices 

0.91+0.08 

1.09 + 0.05 

0.89 ± 0.05 

0.95 

±0.53 

Mid-C, all-hour prices 

1.17 + 0.19 

1.69 ±0.04 

1.95 + 0.04 

1.48 

±0.52 

Mid-C, on-peak prices 

1.21 ±0.11 

1.62 + 0.02 

1.80 + 0.03 

1.40 

±0.52 

Mid-C, off-peak prices 

1.22 + 0.05 

1.33 + 0.04 

1.63 + 0.06 

1.32 

±0.53 


markets according to our preceding study [7]: < a(n) >= 0.90 ± 0.08 for Alberta, < a(n) >- 0.90 ± 0.08 for Mid- 
C. To complement these earlier reported findings, we also conducted a DFA analysis of daily Ontario prices which 
yielded < a(n) >= 0.93 ± 0.08. By comparing these numbers with the < a(n) > values in the table, one can see that 
the discrepancy between the two sets of DFA exponents (the ones obtained for daily and hourly prices) is within the 
statistical uncertainty of our measurements. The fact that the 24-hour aggregation of the prices does not affect the 
average DFA exponent indicates that the latter is insensitive to the short-scale fluctuations of the demand driven by 
the daily socioeconomic dynamics, including the 24-hour cycle and intraday nonstationarities. 

Overall, the correlated price movements in all three markets present an opportunity of forecasting their future 
dynamics based on the historic behaviors over essentially all time scales involved. The Fourier spectral analysis 
confirms the well known existence of daily cycles which can be used as an auxiliary factor in the price forecasting. 

3.2. Probability distributions 

To verify the stability of higher moments of electricity prices, we investigated the Pareto probability distribution 
of hourly, on- and off- electricity data (Fig. 4). The Pareto exponents obtained from this analysis provide a critical 
piece of information on whether the statistical prediction of the prices can yield robust results. 

The probability distributions were estimated using the normalized discrete histograms N(x) with constant bin 
width chosen to be 5 currency units for Alberta and Ontario markets and 1 unit for the Mid-C market. The rate of 
decay of the histogram tails was described by the Pareto index y [47] defined by the power-law fit N(x) ~ x~ (y+1) 
applied to a selected price range. As follows from the definition, the Pareto index describes the asymptotic shape of 
the complementary cumulative distribution function. The values of y were computed using the standard mean-square 
fitting algorithm applied on the log-log scale. 

Depending on the numerical values of the Pareto index, the following cases are possible: if y < 1, neither the av¬ 
erage value of price nor its standard deviation can be determined; if 1 < y < 2, the average price can be calculated but 
the standard deviation cannot; y > 2 means that both parameters can be evaluated and the results of price forecasting 
are potentially robust. The values of the Pareto exponent obtained in our study (Fig. 4, Table 2) indicate that all three 
markets pass the required minimum level y — 2 at least for the price range x e [1,200], Mid-C electricity prices do not 
have a heavy-tailed distribution, so there is no risk of unexpectedly high fluctuations, whereas the Canadian markets, 
especially Alberta, show pronounced fat tails. The probability distributions describing the Alberta market apparently 
has two components - the Poissonian-like central portion with a well-defined maximum, and a heavy power-low tail. 
In Table 2, this behavior is demonstrated by dividing the distributions into two intervals. Pareto indexes of Alberta 
electricity prices in the range x e [200,1000] are lower than in the range x e [1,200] and correspond to the case 
1 < y < 2. Although this does not exclude the possibility of forecasting the average price, such prediction will be 
statistically unstable, especially for prices which are significantly higher than the average. Ontario numbers show the 
opposite tendency, with the range x e [200,1000] exhibiting a faster power low decay compare to the other studied 
price interval. Despite the presence of the fat tail, this market allows for robust statistical prediction even for the 
highest values of the price. 
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Figure 4: Normalized probability histograms of all-hour, on- and off-peak electricity prices. The presented values of the Pareto exponent y indicate 
a possibility of robust statistical prediction of mean price values for all three markets since the condition y > 1 is met. Ontario and Mid-C markets 
are also characterized by stable standard deviations due to y > 2, while the Alberta market does not satisfy this requirement. 


Table 2: Values of the Pareto index y for two ranges of electricity prices in the studied markets. 


Time series 

x e [1, 

200] 

x e [200 

1000] 

Alberta, all-hour prices 

2.072 ± 

0.053 

1.318 

± 

0.068 

Alberta, on-peak prices 

2.360 ± 

0.084 

1.273 

± 

0.070 

Alberta, off-peak prices 

2.394 ± 

0.084 

1.670 

± 

0.082 

Ontario, all-hour prices 

3.016 ± 

0.064 

3.508 

± 

0.045 

Ontario, on-peak prices 

3.177 ± 

0.067 

4.106 

± 

0.128 

Ontario, off-peak prices 

2.964 ± 

0.075 

3.095 

± 

0.115 

Mid-C, all-hour prices 

4.284 ± 

0.190 


- 


Mid-C, on-peak prices 

4.453 ± 

0.155 


- 


Mid-C, off-peak prices 

4.310 + 

0.190 


- 
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All three markets seem to be insensitive to the demand level showing comparable Pareto exponents for all-hour, 
on- and ofF-peak data sets for both price ranges, confirming our earlier observations [7]. It remains to be understood 
why the market demand makes no noticeable contribution to the statistical properties of the process. 


4. Forecasting electricity price movements 


Since anti-persistency implies a stochastic process with negative correlation between its increments, our next 
step is to investigate statistical interdependence of previous and current price movements. To achieve this goal we 
introduced a new method of multiscale increments sensitive to such correlations. The multiscale increments Ax(n),- 
on the scale n are computed according to 



f t~ 1 tj+n > 

J=ti~n -1 t=tj j 


(5) 


and represent a running difference of subsequent aggregated price values obtained by averaging x(t) in non-overlapping 
bins of width n. The statistical relationship Ax, versus Ax, t between the current and preceding increments carries 
information about correlated price dynamics at a given scale and for a given price range. If on average Ax, ~ Ax,_i, 

the price movements are positively correlated and persistent trends are present. On the other hand, if Ax,-Ax,_i, 

the correlation is negative revealing anti-persistent price behavior [29, 7]. 

Fig. 5 shows results of the increment analysis of the three electricity markets for a fixed time scale n — 1. The top 
row of panels displays scatterplots of successive hourly price increments which include both on- and off peak prices. 
The characteristic amplitude of price changes is considerably different for the studied markets, with the increment 
range in Mid-C market being about 10 times smaller than that in Aberta. 

The increment plots of the Canadian markets exhibit a common pattern. Data points in the fourth quadrant of the 
coordinate plane defined by the conditions Ax, < 0 and Ax, [ > 0 are aligned along a diagonal line with slope -1 
suggesting strongly anti-correlated behavior of price increments. 

Mid-C prices also show a correlated regression pattern but with a positive slope revealing positive correlations at 
the studied time scale. Direct comparison of this Mid-C signature with Alberta and Ontario prices is not justified due 
to a drastically smaller range of price fluctuations in the Mid-C market. Because of its compact shape, the regression 
plot of Mid-C electricity prices fits within the core region of Alberta and Ontario markets. 

Besides the main tendency associated with the anti-correlation, Canadian markets show two additional features: 
alignment of the increment values along the positive vertical axis (Ax, > 0, Ax, i » 0) and along the negative 
horizontal axis (Ax, » 0, Ax, | < 0). The first of these features reflects the tendency of the electricity price to 
grow abruptly starting from a steady condition characterized by low price variability. The second feature reflects the 
opposite tendency in which the price drops after a preceding peak and returns to a steady condition. 

The second and third rows of panels in Fig. 5 show scatterplots of price increments constructed for on-peak 
and ofF-peak prices. It can be seen that on- and ofF-peak prices in Alberta and Ontario demonstrate the same anti¬ 
correlation pattern as the one identified for the whole set of hourly prices in these markets. This pattern is more 
pronounced in the on-peak price movements which is an expected result since the amplitude of the price jumps is 
much higher during the peak times. 

Surprisingly, neither on-peak nor ofF-peak price movements in the Mid-C market show positive correlations seen 
in the original all-hour data. A possible explanation of this effect is that the correlations present in the all-hour data 
are produced by transitions between high- and low-demand intervals of the Mid-C market rather than by its dynamics 
during these intervals. 

The last row in Fig. 5 shows aggregated scatterplots constructed by averaging Ax, price increments over a set of 
uniform Ax,_i bins. These are constructed to reveal the prevailing statistical tendencies controlling the price move¬ 
ments in the studied data, enabling their visual and quantitative comparisons. 

It can be seen that the Canadian markets are dominated by the anti-correlated price movements manifested in 
the fourth quadrant of the aggregated scatterplots. The plots exhibit a nearly linear decay, with the slope coefficient 
of about -0.5 for Alberta and -1 for Ontario market. This tendency is observed across the entire range of price 
increments for both on- or ofF-peak data sets. The linear dependence of small price movements was obstructed in the 
original scatterplots but is evident in the aggregated plots. 
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Figure 5: Regression diagrams showing current (A*,) versus preceding (Ajc/_i ) hourly price increments for all-hour, on- and off- peak data. 
Regression plots for Alberta and Ontario markets are similar in shape and reveal negative correlations of price movements associated with the spiky 
structure of the data. The Mid-C price increments exhibit positive correlations in all-hour data only. Bottom panels display averaged versions of 
the diagrams, with Axj values averaged over Ax f _i bins of the same width. 


10 





































The averaged scatterplots of hourly price increments in the Mid-C market suggest positively correlated price 
movements described by a positive plot slop, which is consistent with the shape of the original increment scatterplots 
of this market. To some extent, this tendency is present in the on-peak prices, but not in the off-peak prices. The range 
of price motions underlying this tendency is significantly narrower compared to that in Alberta and Ontario markets. 

The presence of reproducible patterns in the electricity price movements is in a conceptual agreement with our 
DFA analysis results indicating a non-efficient behavior of the studied data described by a + 1.5. Such DFA behavior 
implies that the past and the future price movements are not statistically independent leading to a non-Markovian 
dynamics [47] which could be predicted. However, DFA analysis alone does not provide clues on how to build 
an efficient forecasting model. As we argue below, the specific shape of the increment scatterplots can serve this 
practically important goal. 

It is of particular interest to investigate the uncovered statistical tendencies on different levels of temporal aggre¬ 
gation. Fig. 6 shows scatterplots of subsequent increments between electricity prices averaged over several different 
time sales. As can be seen from the first column of plots, the anticorrelated signature of Alberta price movements is 
manifested on all levels of data aggregation, from several hours to more than a month. The anti-correlated behavior of 
the Ontario market which is evident on the 12-hour time scale weakens quickly with the increase of the aggregation 
scale but can still be recognized even for n = 168 (week-to-week price movements). In contrast, positive correlations 
describing hourly Mid-C prices are essentially absent on all scales exceeding n = 1; moreover, the sign of the correla¬ 
tion changes to negative at n = 12. On this time scale Mid-C becomes qualitatively similar to the Canadian markets, 
speaking in favor of a possibility of universal price movement scenarios in drastically different electricity markets. 

Statistically, the results shown in Fig. 6 indicate that the timing of the price peaks in Alberta electricity market are 
clustered along the time axis, enabling interdependence of price increments across many aggregation scales. Indeed, if 
the price increases occurred at random times, averaging the data would destroy the anti-correlation and the aggregated 
increment plots would show no structure [48]. The clustering effect is much weaker in the Ontario price movements. 
According to our analysis, Mid-C dynamics involves some clustering at the time scale imposed by the daily cycle and 
is not stochastic by its nature. 

The analysis results in Fig. 5-6 confirm that even though the studied electricity markets are deregulated, their 
dynamics can be predicted at many time scales. The forecast in these markets could be built using reproducible 
scenarios shown schematically in Fig. 7. Scenarios I - III represent three characteristic regimes of price dynamics in 
Canadian markets. Scenario IV describes the positive hourly correlations observed in the Mid-C market. It is possible 
that this scenario is also present in the Canadian markets but is masked by their large volatility. 

The scenarios revealed in our study are likely to be applicable to other deregulated markets. Our analysis suggests 
that such markets can function at different levels of price liberalization, with the Mid-C and Alberta markets providing 
the most and the least regulated examples, correspondingly. Ontario represents a stable diversified market representing 
a transitional stage of the liberalization where the floating electricity price does not exceed a reasonable level. The 
asymmetry of the price movements in this market may reflect a specific pricing strategy and decisions, possibly 
working in Canada only. It is widely known that even though the price of the electricity reflects the balance of 
supply and demand, it is also determined by local characteristics of the market such as the availability and the current 
costs of the fuel, the amount of operators, suppliers, and bidders, and other market-specific factors. Therefore, if the 
behavior of the electricity prices is rooted in the management aspects, their regimes can not be universal. The fact 
that both Canadian electricity markets considered here are independent, while the US Mid-C market is actually a pool 
of several submarkets, could be also responsible for the observed differences. It could be the case that the individual 
dynamics of the Mid-C submarkets is more similar to the Canadian markets compared to the entire Mid-C. We leave 
the investigation of this possibility for future research. 
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Figure 6: Diagrams of aggregated price movements for several time scales n ranging from 12 to 720 hours. Most of the plots have a distinctly 
asymmetric shape reflecting a casual relationship between the price movements. For the Alberta and Ontario plots (first and second columns), the 
asymmetry of the cloud of points assumes anti-persistent dependence which can be used for price forecasting. The Mid-C diagrams (third column) 
take different form depending on the aggregation scale, with the persistent and anti-persistent tendencies found at n = 1 (Fig. 5) and n = 12, 
correspondingly. 
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Figure 7: Schematic diagram showing three typical scenarios of price movements. Scenario I: price increase followed by a price decrease (anti- 
persistent behavior); scenario II: stable price followed by price increase; scenario III: price decrease followed by a stable price; scenario IV: 
continuing price increase (persistent behavior). 


5. Concluding remarks 

We have investigated dynamical and statistical properties of price fluctuations in several deregulated electricity 
markets. The results obtained provide new information on the predictability of price movements across a range of 
scale starting from one hour to one month. The following main conclusions have been reached. 

Price fluctuations in deregulated electricity markets are statistically predictable. 

Both methods used for testing the predictability of the electricity prices (scale-dependent DFA and Fourier analy¬ 
ses) clearly show that there are strong anti-persistent correlations between values of electricity prices at all temporal 
scales. This behavior implies that the future and the past price movements are not statistically independent. The time 
series of electricity prices reverse themselves significantly more often than it would be expected for a random walk, 
leading to negative correlations between price increments at all levels of aggregation. 

Statistical forecasts of electricity prices can be statistically stable. 

The obtained values of the Pareto exponent indicate that all three markets are characterized by stable first and 
second statistical moments. This means that the anti-correlated behavior of the markets can be translated into quanti¬ 
tative forecasts with well-defined mathematical expectations and confidence intervals. Even the most volatile Alberta 
market exhibiting heavy-tailed price distributions passes the Pareto test. 

Dynamics of electricity markets can be modeled by several characteristic scenarios working at various time scales. 

The anti-correlated behavior of the studied markets involves a few simple archetypical scenarios which can be used 
for building statistical forecasting models of deregulated electricity prices. The revealed scenarios describe universal 
tendencies of price movements which can be manifested in a different form depending on the time scale of interest 
and the level of market deregulation. 

The presence of negative correlations in the electricity price increments also suggests that they could be simulated 
by a properly chosen mean-reversion stochastic model. The simplest such model can be derived from the Ornstein 
- Uhlenbeck process which combines the standard Brownian motion with the mean-reverting mechanism describing 
the tendency of the price to drift towards a long-term mean supported by fundamentals (see e.g. [47] and references 
therein). More complex mean reversion models, such as the mean reverting jump diffusion model [49, 50], take into 
account sudden price spikes characteristic of deregulated electricity markets. In subsequent studies, the scenarios of 
price movements identified here could be used for adjusting numerical parameters of such models (the mean-reversion 
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level and rate, the volatility of the diffusion component, the frequency of the jumps etc.) as well as for elaborating the 
probabilistic structure of its stochastic constituents in accordance with the non-efficient behavior of electricity prices 
in specific markets. The fact that a does not reach 1.5 at any given time scale suggests that the classical Brownian 
walk approach is not suitable for reproducing the diffusion component of the electricity price movements, and more 
complex stochastic drivers such as those proposed in the fractal market hypothesis [51] would be more realistic. The 
stability of the mean values of the deregulated electricity prices found in our work can be viewed as a new example of 
statistically stable fractal markets [52]. 

In practical terms, these future theoretical efforts would provide an opportunity to reduce commodity and other 
financial risks associated with stochastic dynamics of electricity prices at different time scales using data-derived fore¬ 
casting models. The presented methods provide new relevant information for managing equity derivatives risks over 
time frames exceeding one day. Successful development of this segment of the stock market can make a significant 
contribution to further liberalization process and strengthen the informational efficiency of the electricity market as 
well as the stability of the industry in general. 
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